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Abstract. 

Distant supernovae have been observed to be fainter than what is expected in 
a matter dominated universe. The most likely explanation is that the universe 
is dominated by an energy component with negative pressure ~ dark energy. 
However, there are several astrophysical processes that could, in principle, affect the 
measurements and in order to be able to take advantage of the growing supernova 
statistics, the control of systematic effects is crucial. We discuss two of these; extinction 
due to intergalactic grey dust and dimming due to photon-axion oscillations and show 
how their effect on supernova observations can be constrained using observed quasar 
colours and spectra. For a wide range of intergalactic dust models, we are able to rule 
out any dimming larger than 0.2 magnitudes for a type la supernova at z = 1. The 
corresponding limit for intergalactic Milky Way type dust is 0.03 mag. For the more 
speculative model of photons mixing with axions, we find that the effect is independent 
of photon energy for certain combinations of parameter values and a dimming as large 
as 0.6 magnitudes cannot be ruled out. These effects can have profound implications 
for the possibility of constraining dark energy properties using supernova observations. 



1. Introduction 



Observations of distant type la supernovae (SNIa) show tliat they are fainter than what 
the theory predicts for a universe described by a standard fiat cold dark matter model. 
The most popular explanation is that the universe not only consists of luminous and dark 
matter, but also of large amounts of dark energy, a generic name for a homogeneously 
distributed energy component with negative pressure, p < —p/3\. This is also verified 
by other cosmological probes such as the power spectra obtained from the temperature 
fluctuations in the cosmic microwave background |2] and (on smaller scales) from 
large scale galaxy surveys [3 111 E] ■ 

Several mechanisms have been suggested in order to explain the faintness of distant 
SNIa without having to introduce a dark energy component. Since the existence 
of a large dark energy component has been independently detected by a number of 

I We do not consider any fluctuations of the dark energy component in this paper. 
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cosmological probes, it is unlikely that all the dimming is caused by dust or photon- 
axion oscillations. However, it is crucial to control all possible systematic effects in 
current and future SN surveys focussing on constraining the time evolution of the dark 
energy component. 

Supernova evolution: Some of the apparent dimming of distant SNIa could be 
due to the fact that their luminosities are redshift dependent and that the SNe thus 
cannot be used as standard candles. This could, e.g., be caused by evolution of the SN 
progenitor [01 [Z| • Evolutionary effects can be searched for by comparing the light curve 
shape and spectra of SNe at different redshifts and environments. Most studies show 
that the luminosity properties are rather similar [HI , but others have found variations 
indicating evolutionary effects [HI ITUIITT] . 

Gravitational lensing: Since matter is inhomogeneously distributed in the universe, 
most line of sights to distant objects will probe a matter density significantly lower 
than the average density and the corresponding images will be demagnified relative to 
a homogeneous universe. It is generally accepted that gravitational lensing cannot be 
responsible for all the observed dimming of distant SNIa, however the effect may still lead 
to a non-negligible bias that need to be corrected for, e.g., by using information from the 
magnitude distribution . It has also been suggested that the large scale expansion of 
the universe can be affected by the inhomogeneities in the matter distribution. However, 
perturbative calculations show that this effect is probably negligible [TH| IT^. 

Photon- axion oscillations: It has been proposed that the photon flux from SNIa 
may be diminished because of photons oscillating into axions in much the same way 
as the electron neutrino (z/g) flux, produced in the sun, is diminished because of 
oscillations into muon neutrinos (z/^) and tau neutrinos (z^r)- Due to differences in 
methods and assumptions regarding the magnetic fleld strengths and plasma densities, 
different groups have reached different conclusions about the possible magnitude and 
the energy dependence of the effect [El HEl El UHl 1131111 • 

Intergalactic grey dust: It has been suggested that intergalactic dust could explain 
the observed faintness of high redshift SNIa [211 1221 • there is dust along the line of 
sight, it would extinguish some of the SN light through absorption and scattering. The 
intergalactic dust has to be "grey" in order not to cause a signiflcant reddening that 
would have been discovered if it existed. 

In this paper, we are going to discuss intergalactic dust extinction and photon-axion 
oscillations and show how the effects can be constrained using quasars in the flrst data 
release (DRl) from the Sloan Digital Sky Survey (SDSS). Similar studies have been 
performed using quasars from the SDSS early data release (EDR) for intergalactic dust 
[2Sj and for photon-axion oscillations 24 . For intergalactic dust, a maximum allowed 
amount of dimming of 0.2 magnitudes for SNIa at 2; = 1 was found. In this paper we 
redo the analysis using the larger and higher quality data set of DRl. We also extend 
the range of dust models studied. In the analysis of photon-axion oscillations, there 
was an erroneous factor involving the photon energy in the expression for the density 
matrix evolution. We redo this analysis using the larger data set of DRl with the correct 



Limiting the dimming of distant type la supernovae 



3 



expression for the density matrix evolution. 

In Section |21 we constrain the effects from photon-axion oscillations and in Section 
0] we discuss limits on the intergalactic dust attenuation. In Section |3] we summarize 
our results. 



2. Photon-axion oscillations 



The axion is a hypothetical neutral boson with spin zero originally introduced as part of 
a possible solution to the strong charge-parity (CP) violation. It has also been proposed 
as a candidate for the missing dark matter in the universe. An axion can decay into 
two photons, a — > 27 but also oscillate into a photon. Since they have different spin§, 
photons and axions can only mix in the presence of a mixing agent that preserves the 
quantum numbers such as angular momentum, e.g. an external transverse magnetic or 
electric field 

The photon axion interaction is described by the Lagrangian 

C^-ifE-B. (1) 

where a is the axion field. Ma is a mass scale and E and B are the electric and magnetic 
field, respectively. The mass scale determines the strength of the coupling. Recent 
results from the CERN Axion Solar Telescope (CAST) implies > 8.62 ■ 10^ GeV at 
95% confidence level for an axion mass ma ^ 0.02 eV 



2.1. Theory of the simulations 



To determine the amount of dimming different scenarios of photon-axion oscillations 
lead to, we calculate the conversion probability for oscillations using the density matrix 
formalism. The evolution of a density matrix p, is described by j^ZI 

.dp 



ih 



dt 



(2) 



where 7i is the Hamiltonian, which in our case is given by 7^ = f M. The mixing matrix 
M is given by 123 
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(3) 



where a is the angle in the plane perpendicular to the propagation between a fixed 
polarisation vector and the projected magnetic field. The other quantities in the mixing 
matrix are 
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The axion has spin zero while the photon has spin one. 
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where u is the photon energy, Ue is the electron density, -Bo,± is the strength of the 
magnetic field perpendicular to the propagation of the photon, is the inverse coupling 
between the photon and the axion, and ma is the mass of the axion. The initial condition 
of the equation is 



^1 



Po 



/ i \ 

i 



(5) 

y y 

corresponding to unpolarised light. The diagonal elements correspond to the intensities 
of the two polarisation states of the photon and the intensity of the axion, respectively. 
The system of 9 coupled (complex) differential equations is solved numerically using 
the SuperNova Observation Calculator (SNOC) j2H]- Assuming a quasar (or a SN) at 
a certain redshift, the light beam from the source is followed through a large number 
of spherical cells. Each cell has a value for the magnetic field and the electron density. 
The electron density and the magnetic field strength are evolved as ne oc (1 + z)^ and 
B oc (1 + 2;)^. The dispersion of the field strength and the density is set to 50% of their 
values. We assume that the magnetic field is frozen into the plasma and subsequently 
given by Bq oc nl^^ with random direction ^H]- We also test whether our simulations 
are in accordance with current Faraday rotation measurements. Since the magnetic 
domains are small compared to the total travel length and the magnetic fields have 
random directions, the induced Faraday rotation is too small to be used to rule out 
any of our investigated models. As the beam pass through the cells, the background 
cosmology and the wavelength of the photon are updated, as are the matrices p and M. 
In all of our simulations, we use [^m = 0.3, Q\ = 0.7, h = 0.7]. However, the cosmology 
dependence is weak. 

2.2. Parameter dependence 

The polarisation state of the photon is disregarded by rewriting M as a 2 x 2 matrix, 

where A = A_l = Ay. After solving equation (j2)) analytically we find that the photon 
intensity, pff , is 

P?? = l-§|(l-cosfit), 

1 



n = {A - A^y + 4AI,. (7) 

Using the quantities in equation Q we see that |Am| ^ |A| when ma 3> 
12yT0^ne/(10-8 cm-3)10~^'^ eV ^ 38^nJ{10-^ cm-3)10-^'^ eV. In our parameter 
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interval for Ue, -Bq/Mi, and u, we then have |Am| ^ \^m\ and Q ~ ||^m| oc m\uj ^. 
For ma > 38v/^e/(10-8 

cm~'^)10 eV, oscillations are thus suppressed as m^'^ and 

the dimming is negligible. 

From now on, we assume that ma <^ 38A/^e/(10~'^ cm~^)10~"^^ eV where the 
oscillation probability is independent of the axion mass and the parameters governing 
the oscillation probability are the plasma density, neo, and the magnetic field divided by 
the inverse coupling between the photon and the axion, Bq/M}^^ ||, at 2; = 0. In this case 
we have |A„| < |A| and thus Vt ~ iy/A2 + 4A^. For Buj/{M^ne) > 9 ■ lO^^i cm^G, 
we have |A| ^ 2|Am| and Q ~ Am, i.e., the oscillations are independent of energy. This 
is the most important difference from the results obtained in reference p^, making it 
impossible to exclude any regions to the left of Buj / {Mjif.) ~ 9 ■ 10~^^ cm^G (see, e.g., 
figure IH) . 

2.3. Method 

An effect of photon-axion oscillations is that a dispersion is added to the quasar spectra 
due to the energy dependence of the effect. By comparing the dispersion in observed 
quasar spectra with the dispersion in simulated quasar spectra, we can conclude if the 
model behind each simulation is allowed. 

2.4. Observed quasar spectra 

The spectra of 16 713 quasars are obtained from the SDSS DRl [29j. These quasars 
have redshifts between z = 0.08 and z = 5.41 with a median redshift of z = 1.43. The 
redshift distribution is shown in figure ^ 

The spectra are divided into redshift bins of size Az = 0.2 where we only consider 
the redshift interval 0.1 < 2; < 2.9 since we need a certain number of quasars in each 
bin to obtain acceptable statistics. In this redshift range we have 15 914 quasars. 

Each spectrum is redshifted into its restframe and divided into wavelength bins of 
size 4 A and the average flux is calculated in each of these bins for the individual spectra. 
We calculate the mean spectrum along with the dispersion of the individual spectra from 
the mean spectrum in each redshift bin. All spectra more than three standard deviations 
from the mean are rejected in order to minimize the effect from peculiar spectra. An 
example of a mean spectrum and the spectrum-to-spectrum variation for a particular 
redshift bin is given in figure |21 

2. 5. Simulated quasar spectra 

Using SNOG, we simulate the effects of photon-axion oscillations on quasar observations 
for 10~^^ cm~^ < n^o < 3-10"'^ cm~^ and 10"^^ G < Bq/M^^ < 10"^ G, using the median 
quasar spectrum obtained in reference jHO]- Note that the upper limit of the electron 
density is motivated by the WMAP measurements of the baryon density which gives a 

II In the future we will use the notation M^" = Afa/(10" GeV). 
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Figure 1. Redshift distribution of the quasars in the SDSS DRl. The redshift bins 
in the figure have a width of Az = 0.1. 
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Figure 2. The mean rest frame spectrum and the spectrum-to-spectrum variation for 
the SDSS quasars with redshifts in the interval 1.1 < z < 1.3. 
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Figure 3. Simulated mean quasar spectra for three different scenarios wfien z = I. 

maximum allowed plasma density of Ug = 2.7 ■ 10^^cm^'^|2j. In figureElwe have plotted 
the normalised quasar spectrum for three different oscillation scenarios for z = 1. 

For each redshift bin and set of parameter values, a large number of quasar spectra 
are generated and the mean spectrum and the dispersion due to the oscillations is 
calculated. Comparing the dispersion in the simulated and the observed quasar spectra, 
we are then able to constrain the values of n^o and Bq/M^^. If the simulated dispersion 
is smaller than the observed, we cannot exclude the scenario since real quasars have an 
intrinsic dispersion. If the dispersion is larger for the simulated spectra, a comparison is 
made between the simulated and the observed cumulative distribution of dispersions 
using the Kolmogorov-Smirnov test. The test gives the significance level that the 
distributions differ significantly for the particular redshift bin. The total significance 
level for a certain set of parameter values is calculated as the product of all significance 
levels for the different redshift bins. 

For each scenario we also calculate the rest frame B-band attenuation, due to the 
photon-axion oscillations for a SNIa at z = 1. 

2.6. Result 

In figure |3J it is shown what values of the electron density and the magnetic field 
strength are allowed if there exists an axion with a mass ma < 10~^^ eV. The figure 
also contains the corresponding attenuation for a SNIa. The small scale features of the 
plot is an artefact of the choice of simulated values. If Buo / {MgHe) < 9 ■ 10~^^ cm^G, 
the dimming of SNIa at z = 1 due to oscillations is less than 0.1 magnitudes. If 
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Figure 4. The grey region indicates the allowed parameter region at a confidence 
level of 95%. The lines labeled [0.3, 0.2, 0.1] show the B-band attenuation expressed 
in magnitudes for a SNIa at z = 1. 



Buj / {Mijie) > 9 ■ lO^^^cm^G, the oscillation probability is independent of energy and a 
dimming as large as 0.56 magnitudes for a SNIa dX z = 1 cannot be ruled out. Employing 
the CAST lower limit on Ma [201, the oscillation probability is independent of energy for 
Buj/Uf. > 8- 10~^ cm^eVG. Note that our conclusions are conservative in the respect that 
we have not added any intrinsic dispersion to the simulated quasar spectra. Assuming 
an intrinsic dispersion of 5%, the upper limit on the dimming of SNIa aX z = 1 becomes 
less than 0.05 magnitudes when Bui / [M^rie) < 9 ■ 10~^^ cm^G (see figure EI). 

2. 7. Comparison with the earlier analysis 

The larger sample of the DRl compared to the EDR leads only to minor changes in 
the conclusions. However, the analysis of the effects of photon-axion oscillations using 
EDR data in Mortsell and Goobar |21I contained an error in the evolution of the density 
matrix (corresponding to equation (j21) in this paper). Instead of 2i^ = [7Vl,p], they 
used 2a;i|j = [7Vl,p]. Disregarding the polarisation state of the photons (see Section 
12. 2|) . the photon intensity, p^f , would then be given by^ 

n =^\/(A-A„.)2 + 4AJ,^ (8) 

^ The corresponding equation in reference |24| had also a factor of two wrong due to a misprint. This 
factor was not used in the calculations. 
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Figure 5. The grey region indicates the allowed parameter region at a confidence 
level of 95% when an intrinsic dispersion of 5% has been added to the dispersion 
of the simulated spectra. The thin lines show the B-band attenuation expressed in 
magnitudes for a SNIa at 2: = 1. 

Because of the u in the denominator of the expression for Q, the oscillation probability 
never became independent of wavelength. Therefore, the upper left region of the plot 
could be excluded. 

Another difference compared to the earlier analysis is that = 10~^^ eV was used 
in [21], while ma << 10~^^ eV has been used in this paper. This only affects the region 
where both Ueo and Bq/M^^ are small. 

2.8. Dimming evolution 

If the attenuation caused by photon-axion oscillations is large both at low and high 
redshifts, this cannot explain the faintness of distant SNe, even though it is still possible 
that the oscillations exist and lead to a dimming of all SNe. Thus, we need to investigate 
the axion attenuation as a function of redshift. In figure IHl the difference in observed 
magnitude for different cosmologies is plotted. If the universe has {flu, ^a) = (0.3, 0) 
and the observed SN magnitudes is best fitted by a cosmology with {Qu, ^a) = (0.3, 0.7), 
oscillations need to account for a difference in dimming of 0.20 magnitudes between a 
SNIa at 2 = 0.86 and one a.t z = 0.03. If instead the universe has {flu, ^a) = (1, 0) we 
need to be able to explain a dimming of 0.52 magnitudes between z = 0.86 and z = 0.03. 
In figure in we have plotted the difference in attenuation between z = 0.86 and z = 0.03 
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Figure 6. The difference in magnitude, Am{UM,^A), as a function of redshift for 
different values of SIm and f^A. The sohd Une corresponds to m(0.3, 0.7) — m(0.3, 0) 
and the dotted hne to m(0.3, 0.7) — m{l, 0). 

for different values of tfie parameters rig and Bq/M^^~^. As an example, a universe with 
{Qm,^a) = (0.3,0), He = 5 - 10-^\m^ and Bo/M^^ = 2 ■ lO-^^G where photons mix 
with axions with a mass ma < 10^^^ eV could yield almost exactly the same observed 
magnitudes as a cosmology with (r2M,f^A) = (0.3,0.7) without oscillations (see figure 
ISJ. This scenario is not possible to rule out using the observed (lack of) dispersion in 
quasar spectra. Note however that in this specific case, the oscillation length is smaller 
than the distance to 2; ~ 1. In order to get an attenuation that grows monotonically out 
to 2; ~ 1, we need an oscillation length of the order the Hubble length. For the region 
where the effect is energy independent, this corresponds to Bo/M^^ ~ 10^^^ — lO^^^G. 
We stress again that since cosmic microwave background observations show that the 
universe is close to flat and large scale galaxy surveys indicate a low matter density, 
it is unlikely that all the dimming is caused by photon-axion oscillations. However, as 
for today, it is not possible to rule out substantial systematic effects from photon-axion 
oscillations. 

2.9. Dispersion in type la supernova magnitudes 

Since the oscillation probability is independent of photon energy for Buj/^MsTie) > 
9 ■ 10~^^ cm^G we cannot use the induced dispersion in quasar spectra to constrain 
the effect in this region. As an additional constraint, we can instead use the induced 
dispersion in SNIa magnitudes. If the magnetic field strength and electron density is 

+ Note that the dimming of the simulations are calculated for a universe where (f^M, ^^a) = (0.3, 0.7). 
The effect of the cosmology is, however, small and can be effectively cancelled by small changes in the 
magnetic field strength and the electron density. 
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Figure 7. The solid line indicates the border to the allowed parameter region. The 
different shades indicate regions with seperate differences in attenuation between a 
SNIa at z = 0.86 and one at z = 0.03. The borders to these regions correspond to 
0.15, 0.25 and 0.50 magnitudes. The palest region indicates a difference in attenuation 
of less than 0.15 magnitudes. 




Figure 8. The difference in magnitude, Am, as a function of redshift for {^Im, ^a) = 
(0.3,0.7) (no oscillations) and an open universe with (JImi^^a) — (0.3,0), = 
5 • lO^^^cm'^ and Bq/M^^ = 2 • 10~^*'G. The attenuation is non-monotonic since the 
oscillation length is smaller than the distance to z ~ 1. The error bars correspond to 
the dispersion induced by inhomogeneities in the magnetic field strength and electron 
density (see also Section I^J. 
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Figure 9. The B magnitude and tlie dispersion in the B magnitude as a function of 
redshift for some of the SNe in reference 



not perfectly homogeneous, different lines of sight will correspond to different phases in 
the photon-axion oscillation. This will cause a dispersion also in integrated broad band 
magnitudes of high redshift sources (see figure |H1). Since the dispersion can be of the 
same order as the attenuation, we should be able to constrain the effect by studying the 
observed dispersion in SNIa magnitudes. 

In figure El we plot the B magnitudes of the SNIa in reference [31] as a function 
of redshift. We construct three redshift bins (0 < 2; < 0.1, 0.3 < z < 0.7, and 
0.7 < z < 1.1) with the mean redshifts 0.03, 0.49, and 0.86, respectively. The 
magnitude dispersion in the bins are 0.22 ± 0.02, 0.33 ± 0.04, and 0.30 ± 0.06. We 
use the first redshift bin to estimate the intrinsic dispersion of the SN magnitudes as 
(Tintr = ^/ o"reai ~ "^sim) whcrc (Tfeai IS the dispcrsiou of the observed SNe and cXsim is the 
simulated dispersion caused by photon-axion oscillations at this redshift. We can then 
correct the simulated dispersions in the higher redshift bins for the intrinsic dispersion 
and check if the dispersion is larger than observed. For current data, the dispersion is 
too large to give any useful constraints. If we in the future have a larger SN sample and 
a smaller dispersion we would be able to diminish the allowed parameter region, but we 
would not be able to improve the limits on the dimming, see figure ^1 

3. Dust 

Dust grains are present both in the interplanetary, interstellar, and intracluster space. 
It is not known with certainty if there is dust in the intergalactic space. Dust in the 
Milky Way can be described by a mean extinction law A{X) = f{Rv) where A{X) is the 



Limiting the dimming of distant type la supernovae 



13 




,g-r 

n,o (cm"') 



Figure 10. The solid line indicates the border to the allowed parameter region 
combining current quasar constraints and future constraints using the dispersion in 
SNIa magnitudes. An intrinsic dispersion of 0.1 magnitudes for 50 SNe in each redshift 
bin z = 0.03, 0.49 and 0.86 is assumed. The palest region indicates a difference in 
attenuation of less than 0.15 magnitudes between a SNIa dX z — 0.86 and one at 
z — 0.03. The darker regions indicates a difference of 0.15 — 0.25, 0.25 — 0.50 and 
> 0.50 magnitudes. 



extinction for the wavelength A. Reddening is a consequence of the fact that shorter 
wavelengths are more effectively absorbed and scattered by dust than longer. It is 
quantisised by the reddening parameter, i?y, which is defined as 

^ vy 

where Ay is the V-band extinction and E{B — V) is the colour excess, 

E{B -V)^{B- V)obs - (5 - VU,. (10) 

Here {B — l^)obs is the observed colour and {B — l^)intr is the intrinsic colour. Within the 
Milky Way, we have 2 < Ry < 6, but Ry ~ 3.1 is often a useful approximation. Among 
other galaxies, values of Ry has been observed in the interval 1.5 < Ry < 7.2 |I35| . 
Since small particles preferentially scatter light with short wavelengths, a population 
of small grains have a lower value of Ry than a population of large grains. Different 
parametrisations of f{Ry) is given in, e.g., Cardelli, Clayton and Mathis (CCM) and 
Fitzpatrick [33] • In reference |2S1; the CCM parametrisation was used. In this paper, 
we compare our results using the different forms of f{Ry), finding that the differences 
are negligible except for very high values of Ry. Note however that it is not known if 
any of the parametrisations are valid for an intergalactic dust population with Ry > 7. 
The results presented in this work are obtained using the f{Ry) given in Fitzpatrick 
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3.1. Constraining intergalactic dust 

The main constraints on the intergalactic dust population are: 

(i) The dust population must not produce such an amount of reddening that it causes 
discrepancies with observations. 

(ii) The dust must not contribute with more radiation to the background, than what 
is observed. 

The intergalactic dust absorbs energy from the optical and the UV background and then 
emits it in the far-infrared. This can be used to constrain the intergalactic dust since it 
cannot contribute with more radiation in this wavelength region than what has not been 
identified as other sources. A universe where the dimming of SNIa is explained solely 
by intergalactic dust probably cannot have (f2M7^A) = (I7O), while (i^Mj^A) = (0.2,0) 
cannot be ruled out 

In this paper we will use the reddening of high redshift sources to constrain the 
properties of intergalactic dust and the effects on SN observations (see also reference 
|i36]). This is accomplished with colour comparisons of real quasars with simulated 
quasars for different dust scenarios. 

3.2. Theory of the simulations 

We assume a quasar (or a SN) at a certain redshift and follow the photons emitted from 
the source to the detector at Earth using SNOC |2H!- The path of the photons is divided 
into cells of the same size as the average distance between galaxies. In each cell, the 
dust density and the differential extinction coefficient is specified and the extinction due 
to intergalactic dust is calculated. 

The dust attenuation, Am^^st, at a given emission redshift and a given observed 
wavelength is given by Ao [23] 



where A{V) is the attenuation in the V-band and Dy is the interaction length for photon 
scattering or absorption with dust particles in the V-band. The function H{z) is given 



We study two different models of the dust distribution with pdust oc (1 + z)"^^^ and 
a{z) given by |2Hll2Sj: 

Model A: a{z) = 3 V z 



Model A corresponds to a dust population with a constant comoving density whereas 
model B corresponds to a population where dust is created at redshifts z > 0.5 with the 
same rate as the dilution caused by the expanding universe. The two models represent 




(11) 



by 



H{z) = Hq^{1 + Zf l^M + (1 + Zf (1 - - ^k) + ^^A . 



(12) 
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two extremes in a scale of dust models and they can therefore be considered to capture 
a wide range of possible dust models. 

Since Dv{z) = [crndust(-z)]~^, where a is the interaction cross section and ndust(2;) 
is the dust number density, we have Dv{z) oc pj^g^ oc (1 + z)~°'^^\ The proportionality 
constant is determined by Dov, the interaction length for the V-band at zero redshift. 

3.3. Observed quasar colours 

We use the same sample of 16 713 quasars from the SDSS DRl [29j described in Section 
12.41 The 821 quasars that have been classified as extended are excluded to avoid 
problems with contamination from the host galaxy. 

To avoid quasars with too much noise, we only include objects with magnitudes 
brighter than 22.3, 22.6, 22.7, 22.4, and 20.5, respectively for the five filters (u, g, r, i, 
and z), corresponding to a signal-to-noise (S/N) ratio better than 5:1 ^7\. This leaves 
15 194 quasars. 

We also restrict ourselves to redshifts 0.5 < z < 2. This is primarily because we 
want to have good statistics, i.e., many quasars in every redshift bin. Another reason 
for excluding objects with redshifts z < 0.5 is that quasar spectra at these redshifts 
often are contaminated by their host galaxy. Finally we have 11 694 point-like quasars 
with redshifts 0.5 < ^ < 2 that fulfill the S/N limit. 

The quasars are divided into redshift bins of size Az = 0.05. For each bin, the 
mean colour, {X — Y)°^^, is calculated for all combinations of different filters, X and 
Y. Objects with a colour more than two standard deviations from the mean value is 
rejected. This is because we are not interested in using peculiar quasars that have 
evolved under unusual circumstances. 

3.4. Simulated quasar colours 

We assume a universe with [^m = 0.3,17a = 0.7, h = 0.7] and simulate dust scenarios 
with 10 Gpc < Dqy < 5000 Gpc and < Ry < 12. Each parameter set {Rv,DQy) is 
simulated twice, once for each dust model described in Section 13.21 For each scenario, 
a large number of quasar events are simulated with different values of the redshift. For 
each event, the K correction for the SDSS broadband filters u, g, r, i, and z is generated 
using the median quasar spectrum obtained in reference [2111 • The simulated mean 
colour, (X — yYz'^^^Dov' calculated by subtracting the K corrections for the filters, 
which is the same as subtracting the two filter magnitudes. 

The probability of the different scenarios are determined with the x^-function 

N 

X' = - y)^.,Rv,Dov{V{X - Y)-\,AiX - r).,,K.,Do., (13) 

A(x - r).,«„Do. = (X - r)f ^ - (X - r):x.Dow (i4) 

where i and j refer to the redshift bins. The covariance matrix V{X — Y)ij is defined 
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as 

ViX-Y),, = aiX)l,+aiY)l,+5,, [a(X)?,^ + + a(X - F)^] , (15) 

where crsjs{i) is the systematic error and o"tem(^) is the template error associated with the 
broadband filter i. The template errors arise because we use a median quasar spectrum 
in the simulation of the K corrections and not all rest frame quasar spectra are identical. 
The systematic filter errors are 0.03 for the u and z filters and 0.02 for the others, while 
the template errors have been estimated to 0.1 for the u and z filter and to 0.05 for 
the other filters. (Jobs denotes the deviation from the mean observed colour in each 
redshift bin. Note that the error budget is completely dominated by the template error. 
The inclusion of quasar data from the second and third data release of the SDSS will 
therefore have a small effect on the constraints obtained. 

3.5. Result 

In figure ^3 W6 present our results from the x^-analysis combining A{u — z) and A{r — i) 
together with the expected dimming in the B-band, Amsiz), for a SNIa at the z = 1. 
We conclude that at 99% confidence level, the dimming in the B-band Attlb ^0.1 for 
a SNIa at 2; = 1. We can also limit the interaction length for the V-band at 2; = to 
be larger than 60 Gpc. If the intergalactic dust is similar to the Milky Way dust (with 
Rv ^4), we are able to rule out Am^ > 0.03. 

If we use the COM parametrisation instead of the parametrisation by Fitzpatrick 
we get an upper limit on the B-band dimming of a SNIa at z = 1 of 0.2 magnitudes. 

Note that the scenario without any intergalactic dust is not eliminated by the 
analysis. Also, the increase in the number of quasars and the improvement of the 
systematic filter errors of the SDSS DRl compared to the EDR have led only to minor 
improvements compared to the results obtained in reference |23j . 

4. Summary 

Observations of high redshift SNIa show that the universe is dominated by dark energy 
accelerating the expansion rate. In order to take advantage of the increased SN statistics 
and to be able to use SNIa data to determine the dark energy properties, it is necessary 
to control the systematic effects involved. We have used quasar colours and spectra 
from the SDSS DRl to constrain the effects from intergalactic dust attenuation and 
photon-axion oscillations, respectively. 

We have found that the largest possible dimming of a SNIa at z = 1 caused by 
intergalactic dust is 0.2 magnitudes. The addition of new quasars to the intergalactic 
dust analysis will not improve the limits since the analysis is limited by systematic 
errors, in particular the template errors. One way of improving the template errors is 
to use different template spectra for different groups of quasars. For example, a median 
spectrum created from low redshift quasars will be unaffected by possible intergalactic 
dust and photon-axion oscillations. If the template errors were half their size, we would 
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Figure 11. The confidence levels combining A(u — z) and A(r — i) for the two dust 
models and the attenuation of a SNIa at z = 1. The palest region indicates the 68% 
confidence level allowed by the x^-test and the darker regions indicate 90%, 95%, and 
99%. The almost horizontal lines show the B-band attenuation in magnitudes caused 
by dust for a SNIa at z = 1. 

be able to give a maximum dimming limit of 0.05 magnitudes for a SNIa aX z = 1. 
Also, improving the resolution of the FIRB experiments could lead to the rejection 
of more dust scenarios. The calculations of the dimming caused by intergalactic dust 
is based on mean extinction laws for the dust in the Milky Way [321 EHl- It is not 
known if these parametrisations are valid for an intergalactic dust population where Ry 
might be much larger than in the interstellar medium. Also, we have assumed that the 
reddening parameter, Ry, is independent of direction and redshift. However, this is a 
fair approximation since the simulated SN attenuation is almost independent of Ry (see 
e.g. figure ITTj) . 

For photon-axion oscillations we have found a maximum dimming of 0.1 magnitudes 
for SNIa at z = 1 if Buj / {M^Ue) ^ 9 ■ 10~^^ cm'^G, where B is the magnetic field 
strength, uo is the photon energy. Ma is the inverse coupling between the photon and 
the axion, and is the electron density. However, when Buj / {Marie) > 9 ■ 10^^^ cm'^G 
(corresponding to Buj/Uf. > 8 • 10~^ eVcm^G using the CAST lower limit on Ma |26j). 
the oscillation probability is independent of photon energy and a dimming as large as 0.6 
magnitudes cannot be excluded. We have also shown that for certain combinations of 
and B/Ml^, specifically Bu/iM^rie) > 9 ■ lO'^^ cm^G and Bq/M^^ ~ lO'^^ _ iq-hg, 
the attenuation from oscillations evolve with redshift in a manner very similar to 
the observed SNIa dimming. The region where the effect is independent of photon 
energy could be constrained using the observed dispersion in SNIa magnitudes. An 
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inhomogeneous magnetic field strength and/or electron density will cause different lines 
of sight to correspond to different phases in the photon-axion oscillation. Thus we expect 
a dispersion in integrated broad band magnitudes of high redshift sources if oscillations 
are present. We have shown that it is not possible to use current data to obtain any 
further constraints on this region. However, a future SNIa sample could yield better 
constraints. 
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